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Abstract 

We consider the thermodynamics of the matter made of equal number of neu¬ 
trons and protons and of scalar dibaryons. They interact via the exchange of scalar 
and vector mesons. The interaction is taken into account in the mean field approx¬ 
imation. The condensation of dibaryons in this matter and the phase transition of 
matter to quark matter are considered. Possible signatures of dibaryons in Heavy 
Ion Collisions are speculated on. 
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1 Introduction 


The possible existence and properties of a narrow dibaryon d'(T = 0, J F = 0“) are 
discussed widely in the last years both from experimental [1,2,3] and theoretical [4,5] 
point of view. Being bosons such particles could condense in nuclear matter under proper 
conditions. This phenomenon happens when the nucleon chemical potential grows up to 
half of the mass of the dibaryon. Dibaryons are weakly produced due to the reactions 

n + p —> d! + 7T° 

n + n —» d' + n~ (1.1) 

P + p —» d! + 7T + 

which put them in chemical equilibrium with protons and neutrons. At zero temperature 
the condensation of dibaryons in nuclear (n = p) and neutron matter was studied in 
Refs. [6,7,8]. For the case of d! the condensation starts when the baryon number density p 
grows up to approximately 3p 0 (here p 0 is the normal nuclear density). Such conditions 
(zero temperatures and high densities) are hard to imagine in terrestrial experiments. In 
neutron stars the densities of the central core can be even larger but it is difficult to see 
any clear signal [7] of the dibaryon condensation in neutron stars. 

High densities of nuclear matter can in principle occur in Heavy Ion Collisions (HIC). 
But in that case not only the densities but also the temperatures are high. So we should 
extend the consideration of the dibaryon condensation to finite temperatures. Secondly 
we should take into account the possible phase transition of nuclear matter into quark 
matter. 

Below we will describe our model of nuclear matter with dibaryons in section 2, con¬ 
sider the thermodynamical quantities of nuclear matter in section 3 (condensation curve), 
describe the model for quark matter in section 4 (phase transition) and speculate on 
possible signatures of dibaryons in HIC in the conclusions. 


2 The model of nuclear matter with dibaryons 


We consider an extension of the Walecka model [8] by including dibaryon fields to the 
Lagrangian density 


L = T(iV - m N - g a a - g^u)^ + 


1 2 2 


1 T? 2 I 12 2 , 

4-^V + 2 n V i V+ 


(dt, - - (m D + h a a) 2 ip*ip. 


( 2 . 2 ) 


Here, T is the nucleon field, and a are the c o- and u-meson fields, F ^ — <9 M u j u 

is the field strength tensor of the vector field; tp is the dibaryon field, for which we assume 



that it is a isoscalar-pseudoscalar field. This assumption includes the interesting case of 
the d'-dibar yon. The values m u and m c are the to- and a-meson masses and g w , g a , h 
h a are the coupling constants of the u- and cr-mesons with nucleons (g) and dibaryons 

(h). 

The held operators can be expanded into c-number- and operator parts: 


a = a c + a, 

QfiO^c + Ty, 
P = tp c + if, 
if* = ip* c + f>*. 


(2.3) 


The c-number parts of the fields 71 = <r, uy, <y, and if* are defined as expectation values 
A c =< A > over the ground state of the system. The average values of the operator parts 
are zero by definition: < A >= 0. 

The cr-meson mean held determines the effective nucleon and dibaryon masses in the 
medium 

m* N = m N + g a a c , (2.4) 

m* D = m D + h a a c . (2.5) 

The nucleon vector and scalar densities are defined by expectation values 


Pnv —< v b7oT >, 

Pns =< TT > . 

The vector and scalar density of the dibaryons are dehned by 

Pdv =< P*i do ip- 2h w uj 0 if*if >, 
Pds =< P*P > ■ 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 
(2.9) 


Neglecting the operator parts of the meson fields, we get the following expressions for 
the meson mean fields 

gupNV + huPDV 


r = 


mf 


Or. = 


Pa Pns + h a 2m* D pr>s 


mt 


( 2 . 10 ) 

( 2 . 11 ) 


The nucleon vector and scalar densities are given by 


• y/m*f+k^-u N n 


. \/m*^+k?+v N , 


PNV = 7 J(^(M V,> r^ ) + l)" 1 - (expi ^p- ^ N ) + l)- 1 ), 


( 2 . 12 ) 


Pns = 7 / 


dk 


( 277)3 ^ 

m i?+k" 


((exp ( 7g+gz^ ) + 1} -i + 


(2.13) 



The statistical factor 7 = 4 (2) for nuclear (neutron) matter. Here v n is the nucleon 
chemical potential (excluding electrostatic energy of nucleon in the vector field ) Vn = 

pN ~ 9 ujU C . 

The dibaryon vector and scalar densities are given by 


Pdv = Id I -0p((exp( 


■ V” 


l J+k 2 —Vd 


) — 1 ) 1 — (exp( 


\/m*£+k 2 +v D 


)-l)" 1 ) + 


+Pdv ~ 

= Pdv + P°dvi 


(2.14) 


Pds — Id I 


dh 


( 2jr ) 3 2^/m^+k 2 


((exp( 


- w 1 + (^p( Vm °T +,, ° ) - 1 )-‘)+ 


+ 2^P°DV ~ 
= Pds + P°ds- 


(2.15) 

Here 7 # is the dibaryon statistical factor (1 for d'), Vo is the dibaryon chemical potential 
(excluding electrostatic energy of dibaryon in the vector field) v D = p D — h u (jj C i v D < m * D , 
p° DV is the density of dibaryons in the condensate that is formed when v D = m* D . The 
chemical equilibrium between nucleons and dibaryons means vu = 2 Vn (we consider the 
case when electrostatic energy is not changed in the process of transition of two nucleons 
to dibaryon (h w = 2 g^)). 

The self-consistency equations have the form 


m* N = m N - (g^p ns + K(2m* D p' DS + p° DV )), 


mt 


h, 


m* D = m D -1-( g a pNS + K(2 m* D p' DS + p° DV )). 


m 


(2.16) 

(2.17) 


The total baryon number density equals prv = Pnv + 2 pov- 

We use here the following procedure of solving the self-consistency equations. We start 
with fixing the effective neucleon mass m* N . The value of effective dibaryon mass then 
simply follows from eqs. ( |2.16| , |2T7|) 


m* D = 7 ti d - — ( m N - m* N ). 
9cr 


(2.18) 


Further we assume that the condensation takes place and find the nucleon chemical po¬ 
tential: vn = m* D /2. Then we find the condensate dibaryon density p° DV from eq. (|2.16|) . 
If it is positive the self-consistency equation is solved for a given point. If it is negative 
we take p° DV = 0 and scan v N = v D j 2 from 0 to m* D j 2 to make the eq.( |2.16| ) true. After 
this procedure all the parameters of the matter are defined and we are ready to calculate 
various thermodynamical quantities. 
























3 Thermodynamical quantities of nuclear matter. Con¬ 
densation curve 


The energy density £ =< Too > given by average value of the Too component has the form 


£ = 7 1 (( \J m *N + k 2 + ggj c ) (exp( ^ m * N — ) + 1) 


-l 


+(\J m N + k 2 - g^c)(exp(^ 


l ”f 2+VN ) + l)- 1 ) 


+7d I J0p((yj™§ + k 2 + Ku c )(exp( ^ m * D g k v ~ ) - 1) 

+(\l m D + k 2 - Ku c ){exp( ^ m * D g k —) - l) -1 ) 

+(m* D + Ku c )p° DV + - \mlu 2 c . 


-l 


(3.19) 


The last two terms are here the contributions of the classical u- and cr-meson fields to the 
energy density. 

The hydrostatic pressure p — — | < Ta > has the form 


p = + D- 1 


( 2?r ) 3 V ^/m^+k 2 


i/m^+k 2 

+§7n>/ 

, k 2 , 


[exp( 


k 2 




lv+ k 2 +VN\ -1'1-In 


-) + ir) 


?+k 2 




1 ) 


-1 


V m D + k 2 




1 ) 


-P 


+ 5 m « w c- 


(3.20) 


Because dibaryons in the condensate are at rest they do not contribute to the pressure. 

Now we are in a position to present the results for the thermodynamical quantities of 
nuclear matter. For our purpose the most interesting thing is the condensation curve i.e. 
the curve on the temperature- baryon density plane where the condensation of dibaryons 
starts either at increasing the baryon density or at lowering the temperature. The cal¬ 
culations are made for the standart values of RMF model parameters[9]. The coupling 
constants of the dibaryon <T(2060) to a— and u>— mesons were taken to be[8] 


K 2(yq> h a 1.6g a 


The results are shown on the Fig.l (condensation curve). It means e.g. that at zero 
temperature (T = 0) one has to compress to pxv = 0.54/m~ 3 (~ 2.7 of saturation 
density) to obtain the transition to the dibaryon condensate if the dibaryon mass is 
rri ( ]i = 2060MeI7. Also we show on Fig.l the so called critical curve. It is the curve that 
shows the boundary of applicability of our model. Precisely it is the place where the 
effective mass of the nucleon goes to zero. On the right of the critical curve our model 



cannot be applied. We hope that the more elaborated model will move this curve to the 
right in the unphysical region where the nuclear matter transforms to the quark matter 
for example. And finally we show on Fig.l the phase transition curve i.e. the place where 
the nuclear matter (with or without condensed dibaryons) transforms to quark matter. 
Now we are coming to the description of this curve. 


4 The model of quark matter. Phase transition 


The nuclear matter we are considering is isotopically symmetric (p p = p n ) and has no net 
strangeness (p s = 0 ). At the beginning of the phase transition to quark matter isospin and 
strangeness are conserved. Therefore we will compare the thermodynamical quantities 
of nuclear matter with those of quark matter consisting of equal number of u and d 
quarks and no net strange quarks. The reliable temperatures in our consideration (we are 
interested in nuclear matter with condensed dibaryons) will be below 50 MeV. For such 
temperatures the contribution of pions, kaons etc. to the thermodynamical quantities of 
nuclear matter and the contribution of strange quarks, gluons etc. to the thermodynamical 
quantities of quark matter can be neglected. Then the equation of state (EOS) of nuclear 
matter p N (p N ,T) is described by eq. (|3.20|) and the equation of state of quark matter 
p q (p qi T) is the following one[ll] 


Pg 


' 1 4 i 2rp2\ri 2 ay 16 , ay | 3 


B. 


(4.21) 


ffere ay is the constant of gluonic interaction of quarks. The first term in eq. Q4.21p is the 
pressure of massless u and d quarks, the second term is the plasmon contribution to the 
pressure and the third term is the vaccuum pressure that is responsible for confinement 
(MIT bag model [11]). The constants ay and B were taken to be 


ay = 0.6 


B — 96 MeV/fm 3 


that ensures the stability of massive neutron stars against the transition to strange 
stars[ 12 ]. 

The phase transition of nuclear matter to quark matter starts when the pressure of 
quark matter becomes equal to the pressure of nuclear matter at equal baryon chemical 
potentials and temperatures (Gibbs criterion, here we assume that the phase transition 
is of the first order) 

Pg(p.q, T) p,]y 3 Pq 

The typical dependences of p N and p q on p N for fixed T = 50 MeV are shown on Fig.2. 
The point of crossing of lines is the point of phase transition at given temperature. The 






full phase transition curve is shown on Fig.l. As it follows from Fig.l at temperatures 
above 36 MeV the dibaryon condensation does not occur: increasing the density for fixed 
temperatures above 36 MeV we come to the quark matter rather than to dibaryon con¬ 
densation. The region on the (T, pry) plot where dibaryon condensate is formed (p, n, d! 
in Fig.l) lies therefore below 36MeV for temperatures and above 0.54 1/f m 3 for densities 
(2.7 of normal nuclear density). 

In conclusion let us speculate on the possible signatures of dibaryon condensate for¬ 
mation in HIC. 


5 Conclusion 

The accessible temperatures and densities of nuclear matter in HIC are not very well 
known. Fig.3 shows onr calculations of the temperature of compressed nuclear matter 
formed in HIC as a function of the density of compressed matter (different energies per 
nucleon of colliding ions). We have used the EOS of nuclear matter with dibaryons 
and a one- dimensional shock wave model[13] (Rankin-Hugoniot-Taub equation). The 
compression curve lies above the region of condensation (in ”p,n” or "Q") shown on 
Fig.l. On the first sight we should conclude that the dibaryon condensate never forms 
in HIC. But the compressed matter being formed evolves further and its evolution is 
governed by two effects: by the cooling due to emission of particles (mainly pions) and 
by the decompression (that also results in the cooling). If the radiation cooling prevails 
over the cooling due to the decompression, the matter could evolve into the region of 
the formation of the dibaryon condensate ( p,n,d' in Fig.l). If so, the significant part of 
nucleons will transform to dibaryons (up to 70%). The dibaryons could survive in the 
process of decompression (T^ rs 0.5 MeV) and then decay as free particles in one of the 
following ways 

d! —> nnn + 

d! —» npn 0 (5.22) 

d! —> ppT:~. 

The pions in these decays have the spectrum with the characteristic kinetic energy about 
40 — 50 MeV, if the dibaryon has a mass about 2060MeE. Then these pions could form an 
excess in the temperature spectrum of all pions produced in HIC. In order to these pions 
could be observed they should not be thermalised in the nuclear matter like the other 
pions (direct ones, pions from different other resonances not considered in this paper). 
The principal possibility of this is due to the relatively long lifetime of d'-dibaryon - it 
could decay after the decay of nuclear matter. 



So, we conclude, there is a region of temperatures and densities where the dibaryon 
condensate forms in nuclear matter. If formed in HIC, dibaryons could give a signature 
in form of 40 — 50 MeV kinetic energy pions -the products of their decays. 
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Figure captions 

Fig.l The condensation, phase transition and critical curves on the temperature(T) - 
baryon number density (ptv) plot. The regions of nuclear matter(n,p), nuclear matter 
with condensed dibaryons(n,p, d!) and quark matter(<5) are marked. The critical curve 
shows the boundary of applicability of our model. 

Fig.2 The example of nuclear and quark matter equations of state. The pressures of 
nuclear matter (p n) and quark matter (p q ) are shown as the functions of baryon chemical 
potential(/r at) at the temperature T = 50 MeV. The crossing point is the point of phase 
transition from nuclear matter to quark matter. 

Fig. 3 The temperature of the shock wave(T) as a function of baryon density (p) in ion 
collisions for mean field model of nuclear matter with dibaryons. 



